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Abstract 

' A new dispersion relation for the partial wave vrvr scattering S matrix is set 

■ up. Using the dispersion relation we generalize the single channel unitarity 
;h . condition, SS^ = 1, to the entire complex s plane, which is equivalent to 

the generalized unitarity condition in quantum mechanics. The pole positions 

■ of the a resonance and the /o(980) resonance are estimated based on the 
, theoretical relations we obtained. The central value of the a pole position is 

^ I Mo- ~ 410MeV, To- ~ 550MeV, obtained after including the the constraint of 

^SJ ' the Adler zero condition. 
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PACS numbers: 11.55.Bq, 14.40.Cs, 12.39.Fe 



CN ' Key words: dispersion relation, S matrix theory, chiral perturbation theory 



In a series of recent publications P], 0, H], the present authors presented a dis- 



' persive approach to discuss the single channel and coupled channel tttt interaction 

I I physics. The essence of the method is to use dispersion relations for physical quan- 

titles containing poles and also cuts except those endowed by unitarity. Different 
^ . from more traditional methods, like the K matrix method, Pade approximation, 

etc., in the present scheme the role of the dynamical cuts can be traced explicitly. 
I In Ref. [|], we have established a dispersion relation for sin(25^), where 6t, is the 

' TTTT scattering phase shift defined in the single channel unitarity region. 



sin(2(50 = pF , 

p^(g) = a + V - V - 



^Wlm^^^.^W WO^^.^ (1) 
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where F = ^{S — 1/ S) is the analytic continuation to the real part of the scattering 
T matrix defined in the physical region (times a factor of 2). In Eq. Sj denote 
the possible bound state pole positions and are the corresponding residues of S; 
zY denote the possible resonance pole positions on the second sheet. The integrals 
denote the cut contributions to sin(2(5.n-), L = (— oo,0] is the left hand cut (LHC) 
and R starts from the KK threshold once the 47r cut are neglected, a denotes the 
subtraction constant. The sum of contributions from the cut integrals and from the 
subtraction constant is sometimes called the background contribution. By evaluating 
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the left hand integral in the above dispersion relation using the amplitude (and 

its Pade unitarization) of chiral perturbation theory result, it is found that the LHC 
contribution to sin(257r) is negative and concave in the I=J=0 channel. Therefore 
the o resonance must be introduced to explain the experimental data. 

The present note is a supplementary and an extension to our previous studies. 
It will be shown that the discontinuities of the partial wave S matrix and all other 
physical quantities across the unitarity cut can be expressed as an explicit depen- 
dence on the kinematic factor p — — 'imU s. In other words, the presence of 
the discontinuity on the right is solely due to the presence of the kinematic factor. 
Furthermore wc arc able to re-express the unitarity constraint in a non-trivial and 
analytic expression which holds on the entire complex s plane. For the reason of 
simplicity we will confine our discussion in the single channel unitarity region. As 
an exercise we will estimate the pole positions of the a and the /o(980) mesons using 
the formalism discussed in this note. Some comments related to the coupled channel 
physics will also be made. 

We start from the single channel unitarity region, or more precisely, the center 
of mass energy squared, s, is greater than 4m^ and less than 16m^. The relation 
between the partial wave unitary S matrix and T matrix is defined as 

S{s) = l + 2ip{s)T{s), (2) 
where p — ^1 — Aml./s . With this definition the single channel unitary relation, 

ImT(s) = T(s)p(s)T*(s) , (3) 
is being used together with the property of real analyticity, 

T^*{s + ie) ^T\s-ie) , (4) 

to analytically continue the S matrix and the T matrix, which are analytic functions 
on the physical cut plane, to the second sheet of the Riemann surface: 

T^\s + ze) = T\s - ze) = and S'' ^ . (5) 

One can then verify that the function F defined as 

(6) 

has no discontinuity across the real axis when < s < 16m^, since 

and the left hand cut it contains starts from — oo to 0. The cut structure of F is very 
similar to the cut structure of the function F studied previously. The function F is 
the analytic continuation of cos{25t^) defined in the single channel unitarity region. 
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According to the analytic structure of F, as discussed above we can set up the 
following dispersion relation, 

cos(2(5^) =F = a + J2 



y 2S'{zy){s- z^) tcJl s'-s nJR s'-s ' ^ ' 

where a is the subtraction constant and one subtraction to the cut integrals in the 
above expression is understood. The right hand cut R starts from 16m^ in principle 
but becomes important only when s approaches the KK threshold. [] Using Eqs. (§) 
and (H), we get an analytic expression of S on the complex s plane in terms of poles, 
dynamical cuts, and the kinematic factor: 

S{z) = cos(25^) + i sin(25^) 



2iz-s,) ^ 2pis,)iz - s,) y2p{zY)S'{zy){z-zy) 

+ ^ I ^^ds' + ^ / '^ds' + i / '^ds' + ^ / '^ds'. 

(9) 



TT JL S' — Z IT JL S' — Z TT JR S' — Z TT J R s' — Z 



One may use the definition S'(4m^) = 1 to re-express a in Eqs. @ and (H) in terms 
of other parameters. The above expression respects the well known properties of 
S matrix theory. For example, the physical sheet S{z) does not contain resonance 
poles though the phase motion of S is affected by resonance poles on the second 
sheet. The Eq. (j^), though simple to derive, is an exact relation.^ The Eqs. (||) and 
([^) must satisfy a relation on the whole complex s plane: 

sin^ 25^ + cos^ 25^ = 1, (10) 

which is the analytic continuation of the single channel unitarity relation, S~^S = 1, 
on the complex s plane. The Eq. ([I0| ) is equivalent to the generalized unitarity 
condition S{k)S{k*)* = 1 in quantum mechanics and it contains all information 
about single channel unitarity and analyticity. For example, Eq. must obey 
another relation, 

Siz'/) = 0. (11) 

This equation is derived by the analytic structure of 5* matrix, that is, the physical S 
matrix has zero at the same energy, Zj^ , as the pole energy on the second sheet. The 
Eq. (|11|) is actually equivalent to the requirement of the vanishing of the first order 
pole terms on the l.h.s. of Eq. (|l^) (the second order poles disappear automatically). 
The Eq. (|r^) demands correlations between various parameters: Si, Pi, z^ , S'{zj^), 

^We neglect the An cut from now on in the text. According to the conventional wisdom the 47r 
cut becomes important only above, say, 1.2GeV. 

^Since all cuts at higher energies are actually included in those cut integrals. 
■^We are in debt to the referee who points out this fact to us. 
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a, and cut integrals. However, these relations are in general very complicated to use 
directly. 

We find it helpful, for pedagogical reasons, to analyze Eq. @ together with 
Eq. ( p!OD in some very simple situations. For example, we neglect all the cut integrals 
in Eq. (^), and assume only one pole at s = Sq exist. Then we found two solutions 
satisfying Eqs. and (|TD[): 

1. A bound state: 



a = l-so/2, a = --^so{A- So), /3 = so(4-so). (12) 



The scattering length is a = —y^z^ (taking the mass of the scattering parti- 
cles to be 1). 



2. A virtual state: 
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a = l-so/2, a = -y^so(4 - So), l3 = ^^^^^^ = So(4 - so). (13) 



The scattering length is a 



so 
4-so ' 



We can learn some lessons from these two simple solutions. Comparing with the 
nonrelativistic version of the toy model [Q, 

S = \±^, (14) 

1 — tka 



where k = ys/4 — 1. The nonrelativistic version contains a bound state pole when 
a < 0, and a virtual state pole when a > 0. This agrees with the qualitative behavior 
of the relativistic case. But the pole locates at so = 4(1 — ^) whereas in our case 
the pole locates at So = 4/(1 + ■\). In the nonrelativistic case the pole can locate at 
anywhere between — oo and 4, but in the present case sq G (0, 4). The latter of cousre 
makes sense by eliminating the possible existence of tachyons. Furthermore, in the 
norelativistic case the phase shift 5(oo) goes to ±7r/2 as dictated by the "weak" 
Levinson's theorem. That is not the case in the present situation, even there is no 
dynamical cut.0 All these differences come from the use of relativistic kinematics 
(to use p instead of k), which really makes physical sense, as shown above. The 
relativistic kinematic factor introduces an additional cut from the square root of s 
which is conveniently placed at from to — oo. This additional (kinematical) cut, 
though carefully excluded from the dynamical cuts |1|, 0, does function in its own 
way. For solutions with more than one pole, one can prove that a two-pole (a pair 
of resonances) solution does not exist (in the absence of dynamical cuts), which is 
different from the non-relativistic case. A three-pole solution however exists. As an 
existence proof one can construct the S matrix in the following form: 

s-M^- ipg 

^ = ' 9>0 , (15) 

s — + tpg 



*We call the dynamical cuts as those appear as left hand integrals in Eqs. and (||). 
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which, for > 4 and sufficiently small g, contains a pair of resonances and a 
virtual state pole. In general, however, there is no simple correspondance between 
S matrix poles and the physical resonances |^. 

In the phenomenological discussion on the realistic vrvr scatterings, we follow 
the method of Refs. [|l], to study the properties of resonance poles after the cut 
integrals are estimated. In the following we focus on the IJ=00 channel, using the 
phase shift data from Refs. P, 0, ||, |^. The difference between the fit made in 
Refs. and here is that in Refs. |^ we only fit sin(25^), or the imaginary part of 
the S matrix. In here we ffi the S matrix itself using Eq. (P), 



2i 



ln[S(s)] = Re 



2i 



HSis)] 



+ ilm 



2i 



ln[S(s) 



(16) 



In the single channel unitarity region, 6-^ is real. For the given set of the experimental 
value of 6t,: {sj,6j, A6j}, one may construct the expression of total containing 
two terms, 



xl, in which, 



xl 



\6j - Re 


'iHS{s,)] 


2 







xl 



Im 


" 1 

2i 


HS{s,)] 


2 




A5, 


2 



(17) 



Notice that single channel unitarity in here, unlike most conventional approaches, is 
not guaranteed automatically. If we use the above expression of x^ to make the fit it 
may happen that the x^ minimization program prefers a solution with non-vanishing 
X2, i-e., violating the single channel unitarity. In order to circumvent such a problem 
the unitarity constraint Eq. ( p!OD has to be taken into account to confine the violation 
of unitarity in a numerically acceptable range, which substantially complicates the 
ffi. What we gain with such a price paid is that we can, at least in principle, clearly 
keep track of all kinds of dynamical singularities in their right places. This property 
is not easy to maintain in other approaches which automatically guarantee unitarity. 

The circumvent is possible, noticing that the term x^ in Eq. (0) is in fact quite 
arbitrary, since there is no experimental error bar for the 'imaginary part of 5^'. 
Therefore we can freely chose, for example, another expression of xi; 



^2 ^2 



TlElI"^ 



with sufficiently small e parameter which will guarantee Eq. (0) in a numerically 
satisfiable range. Actually what we do here is an example of the so called 'penalty 
function method' in the theory of probability and statistics . In here the term x^ 
defined in Eq. (0) is called the penalty term and 1/e^ is called the penalty factor. 

In order to make use of Eq. (^ to study the properties of resonance poles, it is 
necessary at ffist to estimate various cut integrals. The discontinuities of function 
F and F on the left can be rewritten as. 



MlF = 2ImLReRT(s) , 
ImL-F = -2p(s)ImLlmRr(s) 



(19) 
(20) 
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since F = 2IiejiT and F = 1 — 2plmjiT. The r.h.s. of Eq. (^) has been esti- 
mated in Ref. [Q, that is one expands ImLReRT(s) to 0{p^) in chiral perturbation 
theory (xPT). But it is easy to see that ImLlmRT(s) vanishes up to 0{p'^) since 

ImRT'^'')(s) = p [ 32^^2 j • Therefore ImLlmRT(s) must be expanded to O(p^), 



ImLlmRT(s) = 2pT(2)lmLRefiT(^)(s) , (21) 

to get a no n- vanishing result. Hence the bad high energy behavior of the chiral 
amplitude gets even worse when estimating Eq. (^), which means when estimating 
the cut-off version of the dispersion integral the numerical result will be very 
sensitive to the cut off parameter. In our understanding the vanishing of ImLlm^T 
at 0{p'^) implies that the quantity and its integral are indeed very small, at least at 
moderately low energies. This suggestion is confirmed by the prediction of the [1,1] 
Fade amplitude which is very small in magnitude.]^ We therefore in the following 
fix the left hand integral of F by using the result from the Fade amplitude. We use 
the same strategy as in Ref. [0] to estimate ImiF and its integral. That is we use 
both the 0{p^) xPT and the Fade approximant to estimate Im^F, for the former we 
truncate the left hand integral at certain scale — which varies within a reasonable 
range. 

One of the lessons one may draw from Ref. p| is that it is not absolutely necessary 
to go to coupled channel situation when discussing, at qualitative level, the property 
of the narrow /o(980) resonance on the second sheet. Therefore we include the /o 
pole in our discussion within the current formalism which only makes use of the 
data in the single channel unitarity region. In some sense, introducing the /o pole in 
the fit improves the determination on the pole location of the a resonance as done 
in Ref. |l[], since in here we no longer need to truncate the data (at around ^/s ~ 
QOOMeV) which is somewhat arbitrary. The right hand cut integrals induced by the 
KK threshold have to be taken into account in here since they will develop a cusp 
structure below the KK threshold. Here we follow the same strategy as in Ref. p| 
to estimate the right hand integrals by using the T matrix parameterization above 



the KK threshold given in Refs. [^Q] and |jTT|, and cut the integral at ^/s ~ 1.5GeV. 
For the integrand we have, 

Im^F = l(ry-l)sin(25^), 
2 rj 

ImnF = ^(--r])cosi25^). (22) 

2p T] 

In fig. [I|, we can see the two estimates of the right hand integral in the IJ=00 
channel. Even though they are not coincide with each other, both of them give the 
same trend when approaching 4m|^. When we only fit sm{2S-,r), the r.h.c. barely 

'^Unlike the situation in the IJ=20 case where the spurious physical sheet resonance (SPSR) 
contribution to cos(2(57r) is large in the IJ=00 case the SPSR contribution to cos(2(5^) is rather 
small. The smallness of the SPSR contribution may be considered a necessary condition for the 
predictions of Pade amplitudes to be numerically reasonable. 
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Figure 1: The contributions from right hand integrals of F and F. Line A is 
obtained from Ref. |T^, hne B is obtained from Ref. |jTl[. 



have any effect to the pole position of /o(980). But in here, we see that the effects 
of the right hand integrals are no longer negligible. 

As already stated earlier the unitarity constraint, Eq.(0), has to be taken into 
account in our fit. Instead of trying to solve the constraints among parameters pro- 
vided by Eq.(|IUp explicitly we make use of the so called 'penalty function' method 
in data fit with constraints among parameters. In principle, increasing the penalty 
factor will drive the fit result moving towards a solution respecting unitarity ex- 
actly. But since there are uncertainties in the input, i.e., the cut contributions and 
the number of pole terms, increasing the penalty factor does not always lead to 
reasonable results. For example, in the present case, for a too large penalty factor 
corresponding to e ~ 0.01 the quality of the fit near the KK threshold becomes 
very bad.[| The masses and widths of the a and the /o(980) poles can be estimated 
from the fit by varying the left and right cut contributions. The variation range of 
the cutoff parameter A in evaluating the left hand integral for sin(25.n-) is taken from 
GOOMeV to SOOMeV here, as we find that larger values of A also lead the fit quality 
below the KK threshold to be rather bad. The e parameter is therefore taken to be 
around 0.02. The results are listed in the following: 

M„ ~ 440 - 530MeV , F^ ~ 540 - 590MeV ; 

~ 976 - 987MeV , Fj^ ^ 22 - 44MeV ; 
a[J ~ 0.230 - 0.276 . (23) 

The above results are compatible with the results of Ref. (the table 1 there), and 
especially Ref. |0] (the Eq. (44) there) though the methods are somewhat different. 
The uncertainty for the width of /o is larger here when comparing with that of 
Ref. , which may be partly due to the fact that in here we only work in the single 
channel unitarity region. In Ref. fl]], the unitarity constraint is not considered since 
only the imaginary part of the S matrix (or sin(2(5^) ) is fitted there. When sin(25^) 
approaches 1, its error behaves as 2co^{25t^)A5t^ and hence approaches 0. Therefore 
the violation of unitarity is automatically confined in an acceptable range in Ref. [|I| . 

^See fig. H, if one takes e = 0.01 the fit curve of (5g would simply miss the data point which is 
just below the KK threshold. 
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From the results we find that the global fit favors a larger value of Oq comparing 
with the results of Refs. @, . A typical fit result is plotted in fig. ||a. The problem 
of having a larger scattering length is due to that we have not put the constraint of 
the Adler zero condition in our data fit. In Ref. this problem is solved by putting 
the constraint of the scattering length parameter by hand. Here we improve the 



(a) 
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(b) 
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g„2 g„2 

Figure 2: (a): A fit using Eq. (^) and the penalty function method. The r.h.c. 
contributions are estimated from Ref. [|T^, the l.h.c integrals are estimated from the 
[1,1] Fade amplitude and e = 0.02. (b): The same as Fig. |^a but with the constraint 
of Adler zero at s = m\/2. 



fit by including the constraint of the Adler zero condition. If we fix the Adler zero 
position at s = wi^/2 the results are the following: 

M„ ^ 380 - 440MeV , F^ ^ 510 - 580MeV ; 

~ 976 - 983MeV , F/„ ~ 43 - 64MeV ; 
a[5 - 0.190 - 0.212 . (24) 

The above results are obtained using the same range of parameters as in obtain- 
ing Eq. (^). The most important change of Eq. (^) comparing with Eq. (|23|) is 
that now the scattering length parameter decreases and is in better agreement with 
the results of Ref. |P, |l3l- We observe again that the decreasing of the scattering 



length parameter drives the a pole moving towards left on the complex s plane, in 
agreement with the observation made in Ref. |]l[. Another major difference between 
Eq. ( PB| ) and is that the latter gives a larger /q width. But Fig. |^ reveals that 



the /o pole is not fitted very well in the latter case. A more reliable determination 
of the /o(980) resonance requires a coupled channel analysis. 

In above discussions one of the major uncertainty in obtaining our results comes 
from the estimates on the left hand cuts which are very difficult to determine ac- 
curately from pure theoretical calculations. Our estimates on the l.h.c. effects are 
based on O(p^) chiral perturbation theory. It means that these estimates are no 
longer trustworthy at high energies, or more precisely, at s ~ and above. In- 
deed, in our calculation to estimate ImF using O(p^) chiral perturbation theory, the 
resonance effects are not taken into account because the resonance only contributes 
at 0{j)^) the real part of the amplitudes. Taking the p resonance as an example, the 
t channel p exchange will enhance the t channel vrvr cut at around t = and will 
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influence the l.h.c. of s cliannel partial wave amplitude through Eq. (53) of Ref. [1]. 
This enhancement will result in a contribution to the l.h.c. starting from around 
4m^ — to further left. If we approximate the effect by a tree level p exchange 
then the effect is approximated by an effective l.h.c. from 4m^ — to — cxd.Q No 
systematic method is known to exist to estimate these resonance contributions to 
the dynamical cuts, in the non-perturbative scheme.^ However, it is reasonable to 
expect that the power counting rule for resonances in perturbation theory [|l^ also 
works in the non-perturbation scheme, at low energies. Therefore it is a reasonable 
speculation that resonance contribution to the left hand integral is 0{p^) at low 
energies and therefore does not distort any of our qualitative conclusion on the cut 
integral at low energy, though their effects will become important at around s = 
and above. Our speculation seems to be supported by a phenomenological analysis 
of Ref. [1^ (see fig. 2 in that paper). A thorough investigation to this problem needs 
a complete new study and goes beyond the scope of the present note. 

To conclude we in this note further extend the previous method we proposed to 
study the partial wave scattering problem by establishing a dispersion relation for 
cos(257r). In our procedure the effects of the unitarity cut are fully exposed by the 
explicit dependence of physical quantities on the kinematic factor. The constraint of 
single channel unitarity, Eq. (|]), is re-expressed as an analytic relation which holds 
on the whole s plane, i.e., Eq. (plQl). The Eq. ([T0| ) is equivalent to the generalized 
unitarity condition in quantum mechanics. Applications of our approach are made 
to determine the pole positions of the a and /o(980) resonances, after estimating 
various cut contributions from both chiral perturbation theory and experiments. 
We find that the central value of the a pole position is about M„ ~ 410MeV, 
To- ~ 550MeV, according to our fit with the constraint of the Adler zero condition. 

One of the author (H.Z.) would like to thank Muneyuki Ishida for helpful discus- 
sions. This work is supported in part by China National Nature Science Foundation 
under grant number 10047003 and 10055003. 
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